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Abstract. We propose to use the threshold-free process of neutrino capture on 
/3-decaying nuclei (NCB) using all available candidate nuclei in the Milky Way as 
target material in order to detect the presence of the Cosmic neutrino background 
(Cz^B). By integrating over the lifetime of the galaxy one might be able to see 
the effect of NCB processes as a slightly eschewed abundance ratio of selected /3- 
^ ' decaying nuclei. First, the candidates must be chosen so that both the mother 

& • and daughter nuclei have a lifetime comparable to that of the Milky Way or the 

signal could be easily washed out by additional decays. Secondly, relic neutrinos 
have so low energy that their de Broglie wavelengths are macroscopic and they 
may therefore scatter coherently on the electronic cloud of the candidate atoms. 
fvj ' One must therefore compare the cross sections for the two processes (induced /3- 

f**'* ■ decay by neutrino capture, and coherent scattering of the neutrinos on atomic 

nuclei) before drawing any conclusions. Finally, the density of target nuclei in the 
galaxy must be calculated. We assume supernovae as the only production source 
and approximate the neutrino density as a homogenous background. Here we 
perform the full calculation for 187 Re and 138 La and find that one needs abundance 
measurements with 24 digit precision in order to detect the effect of relic neutrinos. 
Or alternatively an enhancement of p v by a factor of ~ 10 15 to produce an effect 
p\ ' within the current abundance measurement precision. 



Galactic abundances as a relic neutrino detection scheme 2 

1. Introduction 

In analogy with the photons of the Cosmic Microwave Background (CMB), the 
neutrino component of the early Universe is expected to have decoupled from the hot 
primordial plasma once the rate of interactions with the rest of the particles dropped 
below the Hubble expansion rate [T]. The neutrinos travelling towards us from this 
last scattering surface are normally referred to as the Cosmic Neutrino Background 
(CzaB) or relic neutrinos. 

A detection of the Cz^B would be of great interest to cosmology as well as 
neutrino physics. Contrary to other neutrino sources outside the solar system, 
the Cz/B has the advantage of being a homogenous and isotropic background 
whose number density is only superseded by the CMB photon background. The 
theoretical expected number density is n u+v =112 cm -3 for each neutrino mass state 
[lj. However, this advantage pales somewhat in comparison with the fact that the 
temperature of the background neutrinos is so low that their mean energy is many 
orders of magnitude below the threshold of any of the classical neutrino detection 
methods. 

One can calculate the Cz/B temperature, seeing as it is closely related to the 
very accurately measured photon temperature, T 7 = 2.72548 ± 0.00057 K [2J. Using 
entropy conservation, the photon temperature before reheating - by electron-positron 
annihilation - can be related to the temperature after reheating. The entropy is 
proportional to the number of degrees of freedom times the temperature cubed and 
therefore the final temperature is given by: 

g, tf T} = g.,a*. (1) 

The number of degrees of freedom before reheating comes from both photons and 
electrons, so that g*j = 1 + 7/4. After reheating the only source of radiation is 
photons, i.e. g*j = 1. By assuming the photon temperature before reheating to be 
equal to the neutrino temperature one gets: 

1/3 



Tu = ( - ) T 7 . (2) 



4 
IT 

Consequently the neutrino temperature today corresponds to a neutrino energy of 
only Q.5T^/m u or 3.15T U for non-relativist ic and relativistic neutrinos respectively, 
and a temperature of T v = 1.95K translates to an energy of T u ks ~ 1.6 -10~ 4 eV[3j. 
This energy is far below the threshold for both the radio-chemical and Cherenkov 
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neutrino detectors that have previously been used in experiments with solar, 
atmospheric, reactor, accelerator, and supernova neutrinos. However, it has long 
been known that neutrino capture on /3-decaying nuclei (NCB) does not have a 
threshold (if the normal /3-decay process is already spontaneously occurring in 
nature). Due to energy conservation, the /5-particle will receive all the energy of the 
decay plus one neutrino rest-massji This results in a small peak in the /3-spectrum 
one neutrino mass above the theoretical endpoint, which is basically a smoking gun 
signature for the relic neutrino background unless non-standard model physics is 
introduced. 

The size of the peak is a measure of the number of capture processes, which 
depends of course on the cross section of the reaction, but also on the CjaB number 
density. Recent calculations show that even if one can, in principle, resolve the 
neutrino mass in the /3-spectrum|§|, the background density will have to be enhanced 
by roughly a factor 10 9 in order to produce a detectable CNB signal [4j. Or 
alternatively, the number density of the source material should be enhanced by the 
same factor. 

2. A relic neutrino detection scheme 

Seeing as a O(10 9 ) source-enhancement of any of the upcoming /3-decay experiments 
is not a viable option [5j El [7J [8], we suggest consider the advantages of the NCB 
process in a rather different detection scheme. 

Our approach to the problem consists in letting the detector material for the 
relic neutrino density be basically all the matter in the galaxy. The idea, in short, is 
to look at the effect of neutrino capture on selected nuclei, hoping that the relative 
abundances of selected /3-unstable nuclei - for which both the mother and daughter 
nuclei have extremely long lifetime - will be detectably changed by the presence of 
the neutrino background^]]. That is, the 'detection' would consist in demonstrating an 
anomalous mother- to daughter- nuclei ratio compared with the case without neutrino 

I Under the reasonable assumption that the daughter nucleus is so heavy that its recoil energy is 

negligible. 

§ The calculation in question was performed for the KArlsruhe TRItium Neutrino (KATRIN) 

experiment, which is expected to be the first in a new generation of /3-decay experiments able to 

determine the neutrino mass with sub-eV sensitivity [5]. 

|| This is somewhat reminiscent of the detection scheme presented for supernova neutrinos by 

Lazauskas et al. in |9J 
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induced decays. Put differently, in light of the enormous difficulty of detecting the 
relic neutrino background, we imagine the largest possible detector and calculate 
whether or not this extreme approach will help us achieve the goal. 

The following sections will deal with the details of the calculation. 

Section [3] presents our two proposed candidate nuclei and a calculation of the 
cross sections for NCB on these nuclei can be found in subsection 13.11 However, 
since relic neutrinos have very low momentum and hence a macroscopic de Broglie 
wavelength, coherent scattering on the atomic electron clouds should be taken 
into account. In subsection 13.21 we calculate the relevant cross section for this 
scattering and compare it with the NCB cross section in order to evaluate the relative 
importance of this effect. 

Section H] deals with the production sites and rates for the candidate nuclei must 
be considered. The process contains several important steps: 

• Firstly, which stars are able to produce our choice elements? For simplicity we have 
used supernovae as the only source. It is established that some heavy elements 
can also be produced via radiative pressure processes in e.g. Ap stars [ID], but it 
is fair to say that this is a subdominant site for nucleosynthesis. 

• The next question must be: how much of a given element can be produced in 
a supernova? To answer this we used the nucleosythesis simulations of [11 J and 
produced a simple average abundance plus errorbars, assuming a supernova mass 
range of 10 to 25 M Q . 

• Once we know how much source material is produced per supernova, we need 
to know how many supernovae (or rather progenitor stars massive enough to 
produce a core-collapse supernova) have been produced in our galaxy during its 
lifetime? This question can be answered with a reliable function for the galactic 
Star Formation Rate (SFR). The production rate of supernovae can then be 
constructed in a way that includes the matter distribution of the galaxy in order 
to yield the total distribution of target material. 

• Finally the excess abundance of the daughter nucleus can be calculated by 
folding the target material and neutrino distributions with the cross section and 
performing an integration over time and volume. In addition we need to calculate 
unperturbed abundances. These consists of an integration of the supernova 
distribution (each producing a specified number of both mother and daughter 
nuclei) over volume and time. 
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It should be noted that the CVB is not the only neutrino background in 
our galaxy. The physical background contains also a diffuse supernova neutrino 
background and the diffuse gamma ray burst neutrino background. However, these 
backgrounds are expected to be much more energetic than the relic neutrinos, but 
also much more diffuse. Therefore it is a reasonable first approximation to ignore the 
mixture of neutrino sources and consider the CVB to be the only relevant neutrino 
background. 

3. Nuclear candidates 

For a nuclei to be a suitable candidate in our detection scheme, both the mother 
and daughter nucleus must be very stable, with a lifetime comparable to that of 
the galaxy. Otherwise the signature we are looking for will be diluted by standard 
/3-decays. This requirement in itself presents a problem because the cross section is 
inversely proportional to the product of the Fermi integral and the half life, /ti/2 - 
see Eq. (jSJ) . If we consider the best case scenario of an allowed decay, then log ft 
will typically be in the lower end of the ~ 4 — 9 -interval [13]. Assuming the half 
life to be comparable to the life time of the Universe - i.e. O(10 17 ) s - would then 
imply / to be somewhere in the range of 10~ 8 — 10~ 13 . However, seeing as log / is 
proportional to the Q- value of the /3-decay this leads us to the conclusion that - as a 
second requirement - one must have nuclei with very small decay energies in order to 
achieve both a long half life and a minimal log fti/2 value [T4"l [j~6] . As it turned out 
we were only able to meet the requirement on stability on the mother and daughter 
nuclei. 

In fact we only identified two candidates with both a very long half-life and 
a stable daughter nucleus. These are the first and second order, unique forbidden 
decays of 187 Re and 138 La: 

v e + 187 Re -> 187 0s + e" 
u e + 138 La ->■ 138 Ce + e" 

All additional details are listed in Tables [1] and [2j 
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Decay 


Decay type 


Q- value 


AbundanceM 


Abundance & 


z/+ 187 Re -> 187 0s + e" 


/3- 


2.467 keV 


62.6% 


1.6% 


z/+ 138 La -> 138 Ce + e" 


/3- 


1044.0 keV 


0.0902% 


0.25% 



Table 1: Candidate nuclei for relic neutrino induced beta decay. Data taken from 

d [23 O df . 



Deacy 


log ft 


Tl/2 


Degree of Forbiddenness 


z/+ 187 Re -> 187 0s + e" 


11.28 


1.3727- 10 18 s 


first, unique 


z/+ 138 La -)> 138 Ce + e" 


18.0 


9.3977 -10 18 s 


second, unique 



Table 2: Candidate nuclei for relic neutrino induced beta decay. Data taken from 

/3 03 H3 H3 d- 



5.i. Neutrino-nucleus cross section 

For the calculation of the neutrino-nucleus cross section we have used the following 
expression for the cross section times the neutrino velocity [3]: 

27r 2 ln2 
ctncbVu = —r— — • (3) 

Here A is a description of the ratio of the shape factors involved in the /3-decay 
(denoted 0) and the neutrino capture process (denoted v): 

'm e 



4 



(4) 



C(E e ,p u )„p e E e F(E e ,Z) 

where Wo = m e + Q/3 — rri^ is the maximally available electron energy in the /3-decay, 
F(E e , Z) is the Fermi function and C is the nuclear shape factor. 

In the case of an allowed decay the cross section reduces to the following simple 
expression: 

2tt 2 In 2 



o-ncbVu = p e E e F(E e , Z)'- 



fh/2 



(5) 



And fti/2 can be derived easily from any table of log ft values. 

However, for forbidden decays, the full machinery is needed. According to 
Behrens and Birring, [17| . the shape factor can be simplified somewhat for unique 
K-forbidden processes: 



cm = ^r\ a f£1) 2 u k , 



(6) 
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Neutrino capture on La 
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Figure 1 : The NCB cross section for Rhenium (left) and Lanthanum (right) as calculated 
for three different mass intervals. It is clear that only very large values of m u have any 
visible influence on the cross section of Re, while the La NCB cross section is practically 
not influenced by the neutrino mass at all. Note also that the results presented here are in 
good agreement with Figure 3 of^. 



where the factors preceding uk will be divided out in the expression for A. Meanwhile 
the relevant functions (for our purpose) are given by: 

(7) 



U2 = pt + -^^2p 2 e pl + A 3 Pe- 



The details of the calculations of the coefficients can be found in [17J; for now, suffice 
it to say they depend on the electron energy and momentum and the nuclear charge 
and radius. 

From a computational point of view, the main problem in the calculation of 
Eq. (jBJ) is the presence of r(z)-functions in the A-coefficients, which is to be calculated 
for very small energy input values (that is, very small z). 

The cross sections for neutrino capture on Re and La are presented in Figure Q] 
and agree nicely with the findings of Cocco et al. [3]. 
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3.2. Neutrino- electron cross section 

As mentioned in section [2] we are interested in the cross section for the scattering of 
neutrinos on atomic electrons: 

<s T Abound ' " t & continuum \y) 

To calculate it we use the following expression, due to Bahcall |18| : 

«, = ^/d 3 P | 9 (P)l 2 ^^. (10) 

Here g(p) is the Fourier transform of the bound electron wavefunction and k is the 
four- vector (e e + e u , p e + p„), with (e u , p„) denoting neutrino energy and momentum 
and (e e , p e ) denoting electron energy and momentum. All these quantities are in 
units of m e . Finally the front factor is defined as: 

«b = -(=^r-^ = 1.7xl0-cm'. (11) 



e 



We assume that the scattering reaction of Eq. (Q is only likely to happen 
with the weakest bound electrons. The ground state of atomic 187 Re consists of 
70 electrons in closed orbits and 5 electrons in the 5d (having a binding energy of 
—13.6 eV). We therefore consider the average cross section for scattering with these 
five outer electrons. For 138 La 56 electrons are in closed orbits and the last electron 
is again in the 5d orbital (with a binding energy of -7.3 eV) [19J. 

We now need an expression for the wave function for each of these target 
electrons. In principle a numeric expression can be calculated with e.g. the electronic 
structure modeling programme GAUSSIAN |20j, enabling the user to choose a 
number of specifications, ranging from the class of the basis functions (used when 
expanding the wave function in polynomials) to the stability requirements of say, 
the radius or energy levels in question. However, for easy reference, and because 
we are more interested in the competition between the cross sections and less in the 
precise properties of the atom in question, we have used the Roothan-Hartee-Fock 
wave functions of [19] to calculate the desired Fourier transforms. 

In the Hartree-Fock approach to atomic structure a one-electron wave function 
or spin-orbital, <f), can be constructed as an appropriately weighted sum of Slater 
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function^: 

X = [(2n)\}- 1 / 2 (2C) n+ ^ r"" 1 exp(-Cr)Y lm (0, tp), (12) 

containing an orbital exponent (, quantum numbers n, I, m and the speherical 
coordinates r, d and if. 

As a consistency check we also performed the calculation of the cross section 
using a Thomas- Fermi approximation [21j. We found the two cross sections to agree 
reasonably well, given the nature of the approximations. Figure [2] shows a comparison 
of the calculated cross sections. 

Our results are presented in Figure [3j For Rhenium the NCB cross section is 
fortunately by far the largest for energies less than ~ 10 eV. So we can safely conclude 
that coherent scattering with atomic electrons does not disturb the neutrino capture 
process at the relevant relic neutrino energies. However the second order forbidden 
NCB process on Lanthanum has a cross section of order 10~ 60 which is far below the 
scattering cross section for the entire energy range, and we must conclude that this 
process is simply too weak to be useful in our detection scheme. 

4. Galactic Abundances 

Having calculated the cross sections the NCB interaction rate we now evaluate : 

V(t) = / av v np(z,R,t)dzdR, (13) 



where R and z are the galactic radius and height coordinate, t is the galactic time 
and np(z, R, t) is the number density of the target nucleus. By assuming supernovae 
as the only production sites, the target nuclei density, np, can be expressed as the 
number of target atoms produced per supernova times the supernova rate: 

Tifi = A targotiSN _R SN (14) 

If we consider the neutrino distribution as a homogenous background, our main 
objective - the surplus abundance - can now be calculated with a very simple 
expression: 

^Re+,tot = n u I Tdt. (15) 

^f The total atomic wavefunction is given by the Slater determinant, which can be built from 
a combination of all the electronic wave functions: $ = A((f>\ , ...,<f>n ), where A denotes 
asynnnetrization 
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Figure 2: A comparison of the cross section for Eq. |9]) as calculated with the Roothan 
Hartree Fock method of [19] and a Thomas Fermi approximation for a range of neutrino 
masses. The two results are almost completely on top of each other, with the Thomas 
Fermi result in front. Although the Roothan Hartree Fock cross section is a bit higher than 
the Thomas Fermi cross section the results show a clear overall agreement. The average 
deviation between the two calculations lies between 17.5% (for m v = 0.1 eV) and 16.6% 
(for m u = 1.0 eV). 



The missing ingredient is thus np. In order to calculate it we need firstly an 
expression for the supernova rate. Fukugita and Kawasaki [22] provide the following 
expression: 

J " 0(m)m -1 dm 



R, 



SN 



m 



(16) 



j u <j)(m)dm 
where ip(t) is the star formation rate (SFR) and <f>(m) the initial mass function (IMF). 
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Figure 3: A comparison of the NCB and coherent electron cross sections for Rhenium 
(left) and Lanthanum (right) for a neutrino mass of 0.1 eV. In the case of Rhenium the 
electron scattering becomes subdominant for neutrino energies below ~ 10 eV which is well 
above any realistic relic neutrino energy. However, for Lanthanum the NCB cross section 
is consistently below the scattering cross section and we can conclude that neutrinos are 
far more likely to scatter on the electronic cloud than inducing any second order forbidden 
(3-decays. 



1,35 . As the upper and lower mass limits of Eq. ( fT6l) 
25M , to ensure a mass range that includes core- 



To evaluate i?sN we choose the traditional Salpeter IMF, which for stars heavier 
than the Sun equals cf)(m) ~ m 
we use m c = 1OM and m u = 
collapse supernovae massive enough to actually produce our target materials, yet 
not so massive that the products of nucleosynthesis are swallowed by the ensuing 
collapse into a black hole |23j . With these conservative choices we get: 

R SN = 0.00429M-V0)- (17) 

Following |22] the time dependence of the star formation rate can be expressed as a 
simple exponential law: 

V>(t)=-0(*o)exp[(to-t)/r], (18) 

with r = 2.8 Gyr for z < 1 and t taken to be the age of the Universe. The star 
formation rate is given by: 

log V(*o) [Meyr^Mpc- 3 ] = -2.09+°;^, (19) 
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This expression has been derived in |22| using Eq. (fT6|) and the mean of three 
supernova rates (in units of (10 10 Lbq • lOOyr) -1 ) from three different populations of 
galaxies. To get the local star formation rate we make the following substitutions: 

Wo) -> —~ — J7 Pmw- (20) 

That is, firstly we divide i/j(to) by the local B-band luminosity density of the universe, 
C B = 2.4 ±0.4- lO 8 /iL Mpc -3 , (taken again from [22], where h = 0.72 is the Hubble 
parameter in 100 kms _1 Mpc~ 1 ) and secondly we multiply with the ratio of B-band 
luminosity to stellar mass for the Milky Way (where Mmw,*/Lmw,b = 2.78M /Lb i0 
is taken from [21]). For the matter density distribution we used: 

Pmw(R, z) = — ^ exp(— ). (21) 

*Zd ltd Zd 



Putting these parts together, Eq. (1171) now provides us with an expression that 
describes the physical distribution of supernovae in the Milky Way through a mass 
density function, so that one could in principle combine np with a detailed description 
of the neutrino distribution. 

Specifically our choice of density function is model 8 of Dehnen et al., [25], and 
we have used only the thin disk distribution (which contains most of the matter of 
the galaxy). The following numbers describe the the scale length, Rd, scale height, 
Zd, and central surface density, S^: 

R d = 2400pc, z d = 180pc, S d = 1127.5M p C - 2 (22) 

In total, the star formation and supernova rates are given by: 

i>Mw(R, z, t) = 1.70 x 10- n pMw(i?, z) exp[(£ - t)/r\ (23) 

R SN (R,z,t) = 7.26 xlO" 14 p MW (R,z)exp[(to-t)M. (24) 

Turning now to the second ingredient of Eq. (1T5|) - the amount of target material 
produced per supernova, A^ ta rget,SN - we have used the nucleosynthesis simulations 
of Rauscher et. al. [11] to calculate the average number of target atoms (as well as 
daughter atoms) produced by supernovae in the mass range 10 to 25 M : 

iV Re ,sN = 2.35 ± 0.99 ■ 10 22 iV A 
iV 0s SN = 5.60 ± 2.35 • 10 21 N A 
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Finally we have calculated the unperturbed abundances which is simply the time 
and volume integrated product of /^daughter 1 

A x ,tot= f N x , SN R SN dzdRdt, (25) 

Taking n u = I^+f), Universe to be a completely uniform distribution is not quite 
correct. But as can be seen in e.g. |26] the distribution for a given neutrino mass 
- or analogously momentum - is roughly constant in the inner part of the galaxy, 
which is also where one expects most of the matter to be. With these conventions 
it is very easy to calculate the abundance ratio enhancement in the presence of a 
neutrino density enhancement: 

p _ ^Osncw _ ^Ospid _ ^Os + ^Re+ _ ^Os 
^Re now ^Re i d ^Re — ^Re+ ^Re 

(26) 
We present our final results in Figure HI 

5. Conclusion 

As shown in Figure H] the abundance enhancement is tiny and far below the precision 
of modern galactic abundance measurements even when using meteoritic samples. 
Our results also show that one must have a background enhancement of around 10 15 
to see the effect even on the 8th decimal which is just around the precision available 
on the Osmium abundance - see e.g. |27i [28] 

So despite our use of an extremely large detector volume, we must unfortunately 
conclude that a detection is even more orders of magnitude away with this method, 
than say, with a /3-decay experiment such as KATRIN [5 J or MARE [6j [7J. 
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Figure 4: The Rhenium abundance enhancement as a consequence of the relic neutrino 
background. The left panel shows our standard result, while the right panel demonstrate the 
expected linear enhancement of F^ e as a consequence of a larger neutrino background - in 
the graph p v ^ n indicates the input neutrino density, while p v $ is the theoretically expected 
density corresponding to the normal Uy-y.v—112 cm~^ . Due to the competition with coherent 
scattering processes, we can only assume the semi-constant value of F^ e at low neutrino 
energies to be relevant. Clearly, the standard result of Fr c ~ 1CP 22 is not within reach of 
current abundance measurements and one would need a neutrino background enhancement 
of at least O(10 15 ) in order to actually see the effect of relic neutrinos on the Osmium to 
Rhenium abundance. It can, however, be noted that the error-bars of the plot on the left 
allows for an improvement of approximately one order of magnitude over our mean result. 
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